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ABSTRACT 


We  derive  the  closed  form  expression  for  the  bit  error  probability  of  asynchronous 
dense  WDM  systems  employing  an  external  OOK  modulator.  Our  model  is  based  upon  a 
close  approximation  of  the  optical  Fabry-Perot  filter  in  the  receiver  as  a  single-pole  RC 
fiber  for  signals  that  are  bandlimitr  :  a  ^  equency  band  approximately  equal  to  one 
sixtieth  of  the  Fabry-Perot  filter's  fi*ee  specr  '  ru>ge.  Our  model  can  handle  bit  rates  up  to 
2.5  Gb/s  for  a  fi-ee  spectral  range  of  3800  GHz  and  up  to  5  Gb/s  when  the  power  penalty 


is  1  dB  or  less. 


TABLE  OF  CONTENTS 


L  INTRODUCTION _ 1 

n.  ANALYSIS _ 4 

A.  INPUT  SIGNAL _ 4 

B.  FABRY-PEROT  FILTERED  OUTPUT  SIGNAL. - 6 

C.  DECISION  VARIABLES _ 10 

D.  BIT  ERROR  PROBABILITY _ 15 

HL  NUMERICAL  RESULTS _ 17 

IV.  CONCLUSIONS _ 21 

APPENDIX  A  -  DERIVATION  OF  FORMULA _ 32 

APPENDIX  B  -  PROGRAMS _ 65 

LIST  OF  REFERENCES..... ................. — ..... — ... — .. — ....................... - 70 

INITIAL  DISTRIBUTION. _ 72 


IV 


UST  OF  FIGURES 


Figure  1:  OOK  receiver  structure.................................................... . 22 

Figure  2:  Spectral  characteristics  of  the  Fabry-Perot  filter  and 

the  approximated  single-pole  RC  lowpass  filter.................... . 23 

Figure  3:  Normalized  impube  response  of  the  Fabry-Perot  filter 

and  the  approximated  siuj^e-pole  RC  lowpass  filter................................24 

Figure  4:  Probability  of  bit  error  versus  signal-to-nobe  ratio  as 

a  function  of  normalized  channd  spacing  for  cT  =S...............................25 

Figure  5:  Probability  of  bit  error  versus  signal-to-nobe  ratio  as 

a  function  of  normalized  channd  spacing  for  cT  ■‘10.............................26 

Figure  6:  Power  penalty  versus  normalized  channd  spacing 

as  a  function  of  Fabry-Perot  filter  parameter  cF........... — ....................27 

Figure  7:  Power  penalty  versus  normalized  channd  spacing  as  a  function 

of  Fabry-Perot  filter  parameter  cT  with  a  fixed  threshold  a  =K).5...........28 

Figures:  Worst-case  power  penalty  versus  normalized 

channd  spacing  as  a  function  of  Fabry-Perot  filter 

parameter  cT  with  optimal  threshold  a  =  (A,+A',)/2...............................29 

F^re  9;  Worst-case  power  penalty  versus  normalized 

channd  spacing  as  a  function  of  Fabry-Perot  filter 

parameter  cT  with  fixed  threshold  a  ^  0.5...............................................30 

Figure  10:  The  normalized  optimal  thrmhold  versus  normalized 

channd  spacing  as  a  function  of  Fabiy-Perot  filter  parameter  cr............31 


V 


L  INTRODUCTION 


Asynchronous  wavelength  division  multiplexing  (WDM)  systems  have  been 
increasingly  proposed  as  an  attractive  alternative  to  coherent  optical  frequency  division 
multiplexing  (FDM)  systems  [Ref.  1-S].  Although  asynchronous  WDM  systems  with 
direct  detection  do  not  have  the  channel  capacity  of  coherent  optical  FDM  systems,  they 
are  much  less  costly  to  implement.  Furthermore,  present  filter  technology  enables  the 
designers  to  tightly  pack  the  channel,  resulting  in  asynchronous  dense  WDM  systems  that 
can  provide  aggregate  bit  rates  of  many  terabits  per  second  (1  Tb/s  -  10‘^  b/s). 
Asynchronous  dense  WDM  systems  are  particulariy  attractive  in  the  area  of  undersea 
surveillance  where  hundreds  of  sensors  and  data  collection  sites  are  envisioned  being 
merged  onto  single-fiber  superhighways  through  massive  data  fusion.  Other  applications 
call  for  relatively  low-bandwidth  data  collection  over  many  months  to  be  dumped  quickly 
to  a  remote  recording  she  in  a  matto*  of  minutes.  This  "coUection-and-dumped" 
compression  can  demand  total  data  rates  on  the  order  of  hundreds  of  Gb/s.  Long  distance 
between  data  collection  sites  and  a  remote  recording  site  requires  the  use  of  optical 
amplifiers.  Therefore,  it  is  necessary  to  pack  all  channels  within  the  optical  amplifier 
bandwidth. 

An  asynchronous  dense  WDM  receiver  with  on-ofif-keying  (OOK)  modulation  can 
be  modeled  as  shown  in  Fig.  1.  Conceptually,  the  analysis  involves  two  main  operations; 
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1)  a  convolution  operation  to  evaluate  the  signal  at  the  output  of  the  optical  filter,  a 
Fabry-Perot  (FP)  filter  in  our  investigation,  and  2)  the  integration  of  the  output  of  the 
photodetector.  Evaluation  of  bit  error  probability  by  the  numerical  analysis  of  these  two 
operations  has  been  carried  out  in  [Ref  6],  with  a  number  of  approximations  made  to 
reduce  the  computational  complexity.  In  this  investigation  the  FP  filter  is  shown  to  be 
well  approximated  by  an  RC  filter  within  the  fi'equency  range  l/-/^  |  <  FSR  /  20n,  where 
FSR  is  the  fi’ee  spectral  range  of  the  FP  filter  [Ref  7]  and  /„  is  the  FP  filter  center 
frequency.  For  example,  given  FSR  =  3800  GHz,  the  approximation  works  very  well  for 
\  fifo\'^  effects  of  adjacent  chaimels  within  a  121  GHz  bandwidth 

centered  at  must  be  included,  udiile  all  others  can  be  neglected.  This  simple  model 
agrees  wdl  with  [Ref  6]  as  demonstrated  in  Section  m.  Furthermore,  this  model  enables 
us  to  obtain  a  closed  form  analytical  expresrion  for  the  bit  error  probability  for  which 
numerical  results  can  be  obtained  with  little  effort.  Our  investigation  shows  that  this 
simple  model  provides  accurate  results  as  compared  to  those  in  [Ref  6]  for  bit  rates  up  to 
2.S  Gb/s  when  the  effects  of  two  adjacent  channels  are  included  with  FSR  =  3800  GHz. 
Actually,  whoi  the  power  penalty  relative  to  ringle  channel  operation  is  1  dB  or  less,  there 
is  virtually  no  difference  in  the  effect  of  four  or  two  adjacent  chaimels.  Thus,  for  this 
power  penalty  criterion,  this  simple  model  can  handle  bit  rates  up  to  S  Gb/s  for  a  FP  filter's 
FSR  *3800  GHz. 

In  Section  II  the  closed  form  expression  for  the  decision  variable,  and  consequently, 
the  bit  error  probability  assuming  all  channels  are  bit  asynchronous  as  in  [Ref  6]  is 
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derived.  Section  in  presents  the  numerical  results  which  include  the  bit  error  probability 
versus  the  signal-to-noise  ratio  as  function  of  the  FP  filter  bandwidth  and  channel  spacing, 
and  the  power  penalty  (relative  to  single  channel  operation  without  filtering  or  with 
filtering  but  no  intersymbol  interference)  versus  the  channel  spacing  as  a  function  of  the 
bandwidth.  Finally,  a  summary  of  results  appears  in  Section  IV. 
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n.  ANALYSIS 


The  receiver  model  for  the  asynchronous  dense  WDM  system  is  shown  in  Fig.  1. 
The  desired  signal  is  filtered  by  a  Fabry-Perot  (FP)  filter  that  rejects  adjacent  chaimels. 
The  photodetector  is  assumed  to  have  a  responsivity  (AAV).  The  detected  current  is 
amplified  by  a  low  noise  amplifier  that  contributes  a  postdetection  themud  noise  n{t)  with 
spectral  density  (A^/Hz).  The  decidon  variable  at  the  output  of  the  integration  is 
compared  to  a  threshold  a  to  determine  wh^her  a  bit  zero  or  bit  one  was  present. 

A.  INPUT  SIGNAL 

For  convenience,  we  designate  channel  0  as  the  desired  channel,  and  channel  k  as 
an  adjacent  channel  where  k  =  -M2,  ...,  -1,  1,  ....  M/1  with  M  an  even  integer.  We 
consider  the  equivalent  lowpass  (complex  envdope)  data  signal  in  channel  0  and  channel  k 
as  follows: 

^o(0  =  2  b^jP-iit-iT)  (1) 

/=— Zo 

bk{t)  =  i  ® 

l^—L 

where 

T:  bit  duration 

{0,1}:  bit  in  channel  0  in  the  time  interval  (/T,  (/+1)7) 
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€  { 0,  ^* }  is  the  A*  bit  in  channel  k  in  the  time  interval  (IT  +  t*.  (/  +  1  )r+ 1*) 

;  a  phase  offset  between  channel  k  and  channel  0  and  is  assumed  to  be  uniformly 


distributed  in  (0,  27C)  radians 

(0,, ;  radian  frequency  spacing  between  channel  k  and  channel  0  with  ca„  = 

Xk  :  is  the  time  delay  V>etween  channel  k  and  channel  0  and  is  assumed  to  be 
uniformly  distributed  in  (0,T). 

The  function  - /7)  is  defined  as 

rT\t  ll )  otherwise 

In  both  (1)  and  (2),  the  non-n^ative  int^ers  and  L  represent  the  number  of  bits  in 
channel  0  and  k,  respectively,  that  proceed  the  detected  bits  The  received 
asynchronous  dense  WDM  equivalent  lowpass  signal  at  the  input  of  the  FP  filter  is  given 
by 

r(/)  =  TPbo(fy¥  ^  JPbicif)  (4) 

ic=>^n. 

k*Q 


where  P  is  the  received  optical  power. 
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B.  FABRY-PEROT  FILTERED  OUTPUT  SIGNAL 


The  FP  filter  can  be  chatacterized  by  the  following  equivalent  lowpass  transfer 


function  [Ref  1,7] 


//(/)  = 


t~P 


l-X-p 

1-p 


1-p 


1— j4-~p 
1-p 


(5) 


where  p  is  the  p  -er  reflectivity,  A  is  the  power  absorption  loss  (zero  for  an  ideal  FP 
filter)  and  FSR  is  the  firee  spectral  range.  For  1  / 1  <  FSRJlOn  and  assuming  ^4  =  0,  we  can 
approximate  H  (/)  as  follow; 

H(f)  <=  — = - b- 


=  — b,  tf\<FsmoK 

where 

FSR(l-p) 

C-  p 


(6a) 


(6b) 


The  free  spectral  range  FSR  can  be  related  to  the  full  width  at  half  maximum  (FWHM) 
bandwidth  B  and  the  finesse  F  of  the  FP  filter  as 


FSR  = 


1-p 


=  BF 


(7) 
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Thus  if  the  signal  is  bandlimited  to  |/j  <  FSR/20n,  we  can  truly  approximate  (5)  with  a 
single-pole  RC  filter  with  the  following  transfer  function  and  impulse  response 


//(/)=-W  («) 

1+7-c^ 

h(t)  =  ce~^‘ ,  /  >  0  (9) 

Figures  2a-b  show  the  magnitude  and  phase  (radians)  of  H(f)  of  the  FP  filter  in  (S)  and 
its  single-pole  RC  filter  approximation  given  in  (8)  for  p  =  0.99,  F=3 12.6,  B  =  12. 16 
GHz  and  FSR  =  3800  GHz.  Note  that  as  the  fi’equency  increases,  the  phases  of  the  FP 
filter  and  the  RC  filter  differ  markedly,  but  the  magnitudes  of  their  transfer  functions 
remain  identical  and  attenuate  rapidly.  When  |  /  |  >  FSRJIOk,  the  magnitude  of  H(/)  is 
very  small,  and  therefore,  the  effect  of  adjacent  channel  interference  beyond  this  fi’equency 
range  is  negligible.  Fig.  3  shows  the  normalized  impulse  response  of  both  FP  and 
single-pole  RC  filters.  In  summary,  the  above  approximation  is  valid  for  asynchronous 
dense  WDM  analysis  when  the  filter  finess  F  is  large  or  equivalently  the  FWHM 
bandwidth  B  is  small  since  the  equivalent  lowpass  signal  must  be  bandlimited  to  about 

\f\<Fsmoit. 

This  approximation  has  been  used  in  [Ref  5]  to  study  spectral  efficiency  of  optical 
FDM/ASK  systems,  which  involves  the  evaluation  of  the  decision  variable  for  worst-case 
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analysis  using  the  eye  diagram  technique.  Since  we  are  interested  in  the  detected  bit  b 
in  the  time  interval  (0,7),  we  consider  the  output  filtered  signal  s{t),  0<t<T  given  by 

S{t)  =  SB(t)  +  Sisi(t)  +  SAClit) ,  0<t<T  (10) 


where 

Sg  (t):  desired  signal 

Sjg,(j  ):  intersymbol  interference  (ISI)  signal 
(t):  adjacent  channel  interference  (ACI)  signal 
These  signals  are  evaluated  using  (4)  and  (9)  as  foUows  (detailed  derivations  in  Appendix 
A): 

t 

sb(()  =  /  h{t-X)dK 

0 


=  V?Ao,o(l-e-^0,  0<t<T 


(11) 


=  2  (12) 

/=~Z<o 
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S/icjit) 


=  £  [lit,;  J 

fc=^/2  l=-L 


Jb=-Af/2  /=-t 
h*a 


0</<7’ 


+  6t,-i  / 

-TVCir 


0  <  /  <  X* 


+  i»-i  J 

-T*tt 


Xk<t<T 


+  bt,o  ]  h(t-X)eJ<‘‘<X-^‘^dK] 


Xk<t<T 


=  JP  ‘f  [  ^  it(0. 7)**,;  e-"  -  1) 

h=-MI2  l=-L 


+ii(0,xt)6t,-i  e-^  (e(<H-/«>*y  _ 

’  1+2-r 


+11(T*,  7)6*,-!  e-^ 

c 

+it(Tt,  7)6;to  e'"  (e(«>>*)'  _  (u) 

i'k'J  c 


0<t<T 


where  nitiji)-  {o.  Iilier^e  ^  between  t,  and  tj 


The  FP  filtered  output  s(t)  is  detected  by  the  photodeteaor  which  produces  a 
current  of  fi  1 5(r)  Amps.  This  current  plus  additive  white  postdetection  thermal  noise 
current  fi'om  the  amplifier  is  integrated  by  the  integrator  to  obtain  a  decision  variable  for 
the  threshold  detector. 

C.  DECISION  VARIABLES 

The  decision  variable  Y  appearing  at  the  integrator  output  consists  of  the  signal 
component  X  and  noise  component  N 

Y=X+N  (14) 

where 

T 

X=\  /e\s{t)\^dt  (15) 

0 

T 

N=\n{t)dt  (16) 

0 

We  note  that  is  a  zero  mean  Gaussian  random  variable  with  variance  A'^qT’.  Substituting 

(10)-(13)  into  (IS)  we  obtain  the  signal  component  A* as  a  fimction  of  the  three  parameters 
tjt,cr,  and  a)^^,  which  represent  the  effect  of  intersymbol  interference  and  adjacent 
channel  interference. 

A'=  ^  rnloH  -^(1 

+  ^^(1  -e-^’')[  X 
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+  I  ^  i  btj  i)|2 


lc=-M/2  l=-L  '  l+J-r 

M) 


PT 

+  d’^  Z 

*^M/2 

M/2 

I 

m=-M2 

M) 

ii*M 

pr  M2 

I 

MI2 

2 

m=-M/l 

-I 


©*5^11. 


M2  g(«H/«kX-r«t) 

-MU  m=-Mn  (1+y^Xi-y^)  ^ 


M)  mM 


g(c-/nmX-T*tm)  e(‘^iX-r«k)  g(e-j«>mX-T*tm)  ^ 

+  7^«k  2  J 

*  J  C 


cr  w/o  w«  2(i+y^Xi-;^) 


fc=-M2  I»!=-:M2 

iM  mM 


•>")] 


+ ^  ^  X'  r 

"  ife^-M2  m=^M2 

M)  mi4 


M2 


.r-<Tk  a,k=co„ 


(i+;-Xi-y— )  - 


^  ^  ^*.0  Am.O  gi<^kH  (g-<c^m)T_  g-(Pt>tom)t^) 

k=-M/2  m=-Mf2  ^  *+7^ 


M> 


+ - ;;;5 - i  1 
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+  ^^Re{ 


cT 


T 

k=-MI2 

teO 


M/2 

2 

m=-M/2 


*4,-1  Ko  e(‘^k>"k 


(i+;-rXi-;=^) 


f  g-(<^mytg_g-(o*f»myr  e(c->“*)^m 

L  iTTi"  T  J/ 


1+;=S 


M2  M2  ^  6jy6;;^_, 

^  i  ^  ■“ 


k^-Mn  m^-M/2l=-L 

hM  mM 


(i+y?Ki-y^) 


^  1+y^  2 


cT  I — W«  r  (l+y!^)(l_yl^) 


ifc=^M2  m=-Mai=-L 

hfO  mM 


[ 


1+/ 


bkjb*„g(e~^"*-e~^^ 

_ 


]} 
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r 


+  d’^Re{  T  f 

'fc=-A//2  m=-A//2  (l+;-rXi-y^) 

M>  mM 


_ 

J - e - 


__ 


]} 


(The  above  viiole  term  will  be  zero  if  CD/^  =  Om)- 


i=-Z,o 
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cT  1.  ,  2(1+;^^) 


h^-Ma  l=-L 

M) 


^  Ao,oA*.-i  re'****-!  i-e-<<^*>** 

,  .,.•*  I  .•**  ,  ,•»* 

A— — ^A^72  ^  c  !‘^  c 

M) 


(l_«-^4r)2g(‘^*X-r'«*) 


] 


I  ->  J 

•  •'“  1+y-f  2 


M> 


^  ^0.0^*.0  ^ 
ifc=-M2  1+7?  /? 

M) 


1-7? 


(2e-«*-2e-*^-e-^*4«-2cr)g(«i«t>tt 

2  -I  ^ 


* .  . /^i .  .•***\ 


*=^Ai72  /*-! 

M) 


2(1+7?) 


^<-1  r _  (e«*-e-«*)e-<‘^*>^ 

fc=^W/2  1+7?-  1-7?-  ^ 


M) 


■*  7  J 


+  ^  A*^°  j.£«*-^^*»*e~<‘^*^ 

k^-Mn  1+7? 

M> 


g-^k-e^k  e 


-ZcT 


1-7? 


]} 


(17) 


where  x,  =  min  (x*,  Xm)  “Cg  =  max(xii,  Xm) 
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D.  BIT  ERROR  PROBABILITY 


For  a  detection  threshold  a  and  an  ISI/ACI  bit  pattern  *  =  K  i  =  -L^ 

0;  A=  -A//2,  M2,  0;  the  conditional  bit  error  probal  Uty  of  the  OOK 

signal  represented  by  the  Gaussian  random  variable  Fin  (14)-(16)  is  given  by  [Ref  8] 


where  Q  (X)  is  defined  as 


(18) 


(19) 


and  Jfo  and  X^  are  the  values  of  in  (17)  fi>r  =  0  and  b^ja  =  I,  respectively.  The 
average  bit  error  probability  P,  is  obtained  by  taking  the  expected  value  of  P,  (ft)  in  (18) 
over  all  bit  patterns  ft.  The  minimum  bit  error  probability  is  obtained  by  optimizing  over 
the  thre^old  a.  In  summary,  given  a,  we  calculate  the  following  expectations; 


Pe  =E  {/’,(6)}  (20) 

b 

^ e,  min  =min£  {Peib)} 

O  4 


=  min 


I  jKb) 
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P(6)  =  i(^)"(f)" 


2j*  T 

[/■■■// .-Jgc 

MO  MO 


)d^-Mri--‘d^Mii  dX-Ma  -diwi 


iK 

+j...n 

MO  MO 


Je(^ 


)d^~M/2 .  •  ’d(^M/2  dx-M/i . . .dZMa\ 


(21) 
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in.  NUMERICAL  RESULTS 


In  this  section  we  present  the  numerical  results  for  a)  bit  error  probability  versus 
signal-to-noise  ratio  Z  -  9tP^T/No  as  a  function  of  the  normalized  channel  spacing 
(normalized  to  the  bit  rate) 

=  <“> 

where  4/^=  cOj^/ 2iDb  is  the  equal  channel  spadng  in  Hz,  andb)  power  penalty  versus 

normalized  channel  spacing  7  as  a  function  of  the  FP  filter  parameter  cT - 

where  B  is  the  filter  full  width  at  half  maximum  bandwidth(FWHM).  For  our  model  to  be 

valid  we  set  A/  -  2,  that  is,  we  constraint  the  signal  to  be  bandlimited  to  within  two 

adjacent  channels,  thus  we  incorporate  the  degradation  the  two  nearest  adjacent 

channels.  We  observe  that  for  bit  rates  of  1  Gb/s  or  less,  oui  model  is  valid  up  to  A7  ==  10 

and  veiy  little  difference  is  observed  between  M  =  2  and  M-  10.  Also,  we  observe  that 
there  is  little  difference  betweoi  M=2  and  A/=  4  when  I>  10  for  bit  rates  up  to  3  Gb/s. 
In  all  results  we  set  Lg  =  2  and  L  »  1. 

Figures  4-5  show  the  minimum  bit  error  probability  versus  signal-to-noise  ratio 

(  Z )  for  cT  =  5  and  10,  respectively.  Along  with  each  curve,  we  also  show  that  of  the 
synchronous  case,  i,e,  X*  =  0,  =  2  and  L  =  0,  and  that  of  a  single  channel  (SC) 

operation  without  filtering  or  with  filtering  but  without  ISI.  In  Fig.  4  we  observe  that  a 

large  degradation  occurs  due  to  ISI  for  cT=  5  which  represents  a  narrowband  filter.  As 
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the  FP  filter  bandwidth  is  made  larger  as  in  Fig.  S  with  cT=  10,  the  ISI  is  reduced  but  the 
ACI  increases. 

In  our  model,  we  are  constrained  to  Af  =  2  for  the  case  under  consideration.  We 

use  an  FP  filter  with  FWHM  B=12,16  GHz,  firee  spectral  range  FSR  =  3800  GHz,  finess 
-p)  =  FSRIB  =  312.6,  and  c  *  38.4  GHz;  1/7  *  2.56  Gb/s,  and  cJ  =  15.  If 
the  powa  penalty  related  to  single  channel  operation  is  to  be  less  than  1  dB,  then  by 

Fig  6,  a  minimum  normalized  chaimel  spacing  of  /  =  12  must  be  used.  Equivalently 

from  [Eq.  (22)],  we  could  use  4^  »  7/7  »  12*2.56  GHz  s  30.4  GHz,  (Le.  the  channel 

pacing  is  a  multiple  of  the  bit  rate).  In  this  model  the  farthest  adajcent  channel  for  2 

is  twice  the  channel  spacing  which  is  60.8  GHz.  This  verifies  the  assumption  1  /-/g  1  < 

FSR/IOk  s  60.5  GHz,  where /g  is  the  FP  filter  center  frequency.  This  result  agrees  well 

with  that  in  [Ref.  6:  Figs.  6,9,  M/F  »  0.4,  a  =  0 Ji].  Thus  we  incorporate  the  degradation 

caused  by  the  two  nearest  adjacent  channels.  We  observe  that  for  bit  rates  of  1  Gb/s  or 

less,  our  model  is  valid  up  to  Af  =  10,  and  very  little  difference  is  observed  between  Af  = 

2  and  Af  s  lO.  Also,  we  observe  that  there  is  little  difference  between  Af  «  2  and  Af  »  4 
when  /  ^  10  for  bit  rates  up  to  3  Gb/s.  In  all  results  we  setLg  =  2  and  L  »  1. 

Figure  6  also  shows  the  power  penalty  for  an  asynchronous  dense  WDM  system 

relative  to  a  single  chaimd  operation  at  the  minimum  bit  error  probability  of  10''^  This  is 

the  required  additional  signal  power  (dBW)  for  the  asynchronous  dense  WDM  system  to 

be  able  to  operate  at  the  10''^  bit  error  probability  achieved  in  the  single  channel  system 

with  a  SNR*12dB.  The  asynchronous  dense  WDM  system  is  /57-limited  at  2. 15  dB,  0.9S 

dB,  0.6  dB,  and  0.5  dB  in  power  penalty  for  cT  -  5,  10,  15,  and  20,  respectively.  It  is 
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seen  that  for  a  2.3  dB  power  penalty,  the  normalized  channel  spacing  can  be  as  close  as  / 

=  6  (i.e.,  a  channel  spacing  of  six  times  the  bit  rate)  for  cJ  =  5.  If  the  power  penalty 
criterion  is  1  dB,  the  normalized  channel  spacing  is/»  12  for  c7’=  10,  IS,  20.  We  remark 
that  although  the  exact  transfer  functions  of  the  FP  filter  is  used  in  [Ref  6],  a  number  of 

approximation  have  been  made  to  obtain  numerical  results.  The  approximations  are  1) 

the  ISI  is  obtained  by  modeling  FP  filter  as  a  single-pole  RC  filter  [Ref  6,  Eqs.  (4)  and 

(36)],  2)  approximating  the  finite  integration  with  an  infinite  integration  in  the  calculation 

of  ACI  [Ref  6,  Eq  (IS)],  and  3)  the  beat  interference  is  ignored.  On  the  other  hand,  the 

ISI  and  ACI  in  our  investigation  are  obtained  by  modding  the  FP  filter  as  a  single-pole  RC 

filter,  using  finite  integration  and  including  the  beat  interference.  Since  the  results  in  our 

investigation  and  in  [Ref  6]  agree  well,  we  conclude  that  approximations  are  quite  valid. 

We  also  note  that  our  results  also  t^ee  well  with  the  simulation  carried  out  in  [Ref  1, 

Fig.  17]. 

The  above  numerical  results  shown  in  Figs.  4-6  are  obtained  with  an  optimized 
threshold  setting.  Figure  7  shows  the  power  penalty  for  fixed  threshold  a  =  £PTI2  which 
is  the  same  optimum  threshold  for  single  channel  operation  (midpoint  between  the 
received  power  for  bit  zero  and  bit  one).  It  is  seen  that  the  performance  of  an 
asynchronous  dense  WDM  system  is  quite  sensitive  to  a  for  a  narrow  band  filter.  An 
additional  1.5  dB  is  observed  for  cJ  =  5  for  /  >  8,  and  0.4  dB  for  cT  =  10  for  /  >  12. 
Negligible  degradation  is  observed  for  15, 20  for  I  >  16. 
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Figures  8-9  show  the  power  penalty  versus  normalized  channel  spacing  as  a  function 
of  FP  filter  parameter  cT  for  the  worst-case  analysis  with  optimal  threshold  and  fixed 
threshold,  respectively.  The  worst-case  bit  pattern  is  fixed  to  produce  the  minimum  A', 
and  maximum  Aq  where  A,  andA,,  are  the  values  of  A  in  Appendix- A  equation  (4-21) 
with  ^0.0 “  1  ^0.0 “ 0  respectively. 

We  observe  that  the  power  penalty  for  the  worst-case  analysis  is  only  slightly  larger 
than  that  of  the  exact  analysis  for  /  >  10  shown  in  Fig.  6.  Similarly  the  power  penalty  for 
the  worst-case  analysis  with  fixed  threshold  is  only  slightly  larger  than  that  of  the  exact 
analysis  with  fixed  threshold  for  /  >  10  shown  is  Fig.  7.  The  reason  for  this  is  that  for 
large  channel  spacing  (/  >  10),  the  ACI  effort  is  small,  so  the  ACI  bit  pattern  has  a  small 
influence  on  the  power  penalty. 

Figure  10  shows  fue  normalized  optimal  threshold  for  the  exact  analysis  shown  in 
Fig.  6.  It  is  observed  that  a  »0.4  for  I  >  10.  Note  that  the  normalized  optimal  threshold 
for  the  single  channel  operation  is  a  =  O.S. 
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IV.  CONCLUSIONS 


We  have  presented  a  simple  model  for  the  analysis  of  asynchronous  dense  WDM 

systems  employing  an  external  OOK  modulator.  The  only  approximation  that  we  use 

involves  the  modeling  of  the  Fabry-Perot  filter  by  a  single-pole  RC  filter  assuming  the 
equivalent  lowpass  signal  is  bandlimited  to  the  fi-equency  range  I/I  <FSR/20n.  This 
model  enables  us  to  obtain  a  closed  form  expression  for  the  bit  error  probability  which 

previously  could  only  be  obtained  via  numerical  analysis  [Ref  6].  For  an  FP  filter  with  an 

FSR  around  3800  GHz,  our  model  can  include  the  ACI  effects  of  two  adjacent  channels 

for  bit  rates  up  to  2.5  Gb/s.  Our  numerical  results  show  that  this  model  agrees  well  with 

that  in  [Ref  6]. 
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Figure  I  :OOK  receiver  structure 
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Figure  2:  Spectral  characteiistics  of  the  Fabry-Perot  filter  and  the  ^proximated 
single-pole  RC  lowpass  filter 
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time  (ns) 


Figure  3:  Normalized  rimpulse  response  of  the  Fabry-Perot  filter  and  the  approximated 
single-pole  RC  lowpass  filter 
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Power  Penalty  (dB) 


Figure  6:  Power  poiahy  versus  normalized  channel  spadng  as  a  function  of  Fabiy-Perot 
filter  parameter  cT 
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Figure  7:  Power  penalty  versus  normalized  channel  spacing  as  a  ilinction  of  Fabry-Perot 
filter  parameter  cT  with  a  fixed  threshold  a^^O.S 
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Figure  8:  Power  penalty  normalized  channel  spacing  as  a  function  of  Fabry-Perot 
filter  parameter  cT  for  worst-case  with  optimal  threshold 
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Figure  9:  Power  penalty  versus  normalized  channel  spacing  as  a  function  of  Fabry-Perot 
filter  parameter  cT for  worst-case  vdth  fixed  threshold  a=O.S 
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optimal  Thrsahokl 


Figure  10:  The  normalized  optimal  threshold  versus  normalized  channel  spacing  as  a 
function  of  Fabry-Perot  filter  parameter  cT 
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APPENDIX  A 


Derivation  of  Formula 


/.  INTRODUCTION 


Desired  channel 
Adjacent  channel  : 


Channel  0. 

Channel  k. 

A:  =  -  M/2, ...  -1,  1, ...  M/2,  M :  even. 


Bit  in  channel  0  in  time  interval  (/T,  (i+l)T)  :  e  {  0, 1 } 

Detected  bit  in  channel  0  in  (0,T) :  ^  e  {0, 1 } 

Bit  in  channel  k  'ml^  time  interval  {IT  +  Tit,  Q  +  1)7’  +  X/t)  ' 
bkj  e  {0,e'**},7= 

frequency  spacing  between  channel  k  and  channel  0, 0)jt  =  -(O-jt. 

(|)^  :  phase  offret  between  channel  k  and  channel  0,  uniformly 
distributed  between  (0, 27c). 

Tjt :  time  delay  between  channel  A:  and  channel  0,  0<Tii<T. 


Data  signal  in  channel  0  :  bo{t)  =  2  Prit  —  iT) 

»=-Lo 

0 

Data  signal  in  channel  k :  bk{i)  =  2  Prit  -  (/7’+  Tit)) 

l=-L 

Lq,  L  :  positive  integers. 

P 

-*  Tvv  1  0, otherwise 

The  received  signal  at  the  input  of  the  FP  filto’  of  channel  0  is  : 

r(t)=JPbo(i}+  ^  JPh(t) 

k=-MI2 

M> 

P :  received  optical  power. 


(1-1) 


(1-2) 
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The  V  tput  of  the  FP  filter  is: 


00 

'‘0(0  =  J  /t(/  -  X)r(k)dk 

'-00 


r—  r  r—  m/2  r 

=  Jp\Kt-X)b(,(X)dX+JP  £  \  h(t-X)bi,(X)dX 

-00  k=-MI2  -00 

il[*0 

00 

=  V?6o,oJ  h(t-X)PT(X)dX 

-00 

-1  ® 

+  JP't  bo,i  j  h{t-X)PiiX-tT)dX 

—00 

+/J“  ‘f.  {  i  bn  1  h(t-X)e^^^-'‘^PTQi-(lT+xt))dX 

k=-Ma  h-L  ’  -<o 

fato 

00 

+6*,-,  /  A(/  -  -(-?’+ 

-00 

00 

+6t.oJ  hO-X)e>‘“^^-'‘’^PT(.X-Xk)dX  } 

-00 


h(t)  is  file  equivalent  lowpass  impulse  response  of  the  FP  filter  of  channel  0. 
Smce  the  detection  intend  is  0^<T,  we  only  need  to  evaluate 

=  ^0  (f)>  0<t<T  and  write  S(tJ  as  : 

S(i)  =S/t)+  SJt)  o<t<T 

Sg  (t)  :  Desired  signal  in  channel  0,  the  1st  term  {i=0)  of  (1-3) 
S^(t)  ;  Inter-symbol  interference,  the  2nd  term  (i  <  -1)  of  (1-3) 
Adjacent  channel  interference,  the  3rd  term  (E)  of  (1-3) 

k 


All  evaluated  in  0<t<T. 


//.  CALCULATE 


In  this  section  we  are  to  calculate  the  easy  part  of  S(t)  —  Sg  (t)  and  (t). 


A.  FP  filter : 


The  lowpass  equivalent  transfer  function  of  FP  filter  is 


m- 


l-p  l-A-p  1-p  \-A-p 


p  power  reflectivity. 

A  power  absorption  loss  (0  for  ideal  filter). 

FSR  :  free  spectral  range. 

Since  f«FSR  (for  operating  fi-equency  range),  we  can  approximate 
H(f)  2s  (assume  A  =  0): 


l-p  l-p 

”  /I  _  W  ^  .  .•  2n/p 

(1  PH/PfSR  (l-p)FSR 


l+/-r 


where  c  = 


FSR(l-p) 


for  FP  filter,  we  also  have  • 

rcD  "^Jp 

,  B  :  full  width  at  half  maximum  bandwidth  (FWHM) 
of  half  power  bandwidth. 

Using  inverse  Fourier  transform,  we  find 
"\t)  “10,  otherwise 
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B,  Calculate  ft) : 

00 

S,(t)  =JPbo,o  f  h(t-\)PT{X)dk 

-00 

t 

=  JPbdfi  J  h{t  -X)dk 

0 

t 

=  JPbaj:)  J  dk 

0 

/ 

=  JPbofi  ce^"  j  dk 

0 

=  v'?6o.o(l- O 


(a)  For  Bit  "1"  :  8^,(1)  ~  0<t<T 

For  Bit  ”0”  :  Sgo(t)  =  0  0<t<T 


C  Calculate  S^(t) 


-1 


00 


S^(t)=jp  i:  V  f  h(t-k)PT(X-iT)dK 

i=-Lo 


-1 


-00 
(i+nr 


=  JP  t,  V-  J  ¥t->-)dk 

i=-Lo  iT 

-1 

=  VF  2  ^0,/  J  ijfX 

i=-Lo  jT- 


-1 


=  JPe-^  2  bo^  -  e‘‘^)  dX 

i=-Lo 


(2-1) 


0<t<T 


(2-2) 


IIL  CALCULATE  S^t) 

In  this  section  we  are  to  find  with  the  procedure  ; 

A.  When  /  <  -2 

B.  When  /  =  - 1 

C.  When  /  =  0 

D.  Summary  for  5^^  0) 

A,  When  I  < -2 

Integrate  the  first  term  of  b^(tj  in  (1-3),  we  have  ; 

_  M/2  -2  “2 

Jp  I,  H  bkj  I  Kt  -  X)  -  (.IT+  ^k))dX 

k=-M/2  l=~L  Jqo 

M> 

_2  (/+i)r+Tt 

=  JP'L'tbn  f  h{t-X) 

.2  ('+•)?■•«» 

=  JP'ZI,bkj  f  ce-^'-« 

^  ^  /r-w* 

_2  (/+i)r+tt 

=  y?226*.;  I 

^  ^  ’  /rix* 

g-{ct+j<akTk)g{c+j<Ok)ilT-i^k)^gic+j^  _ 

it  /  1+/^ 

=  JP  1.1,  bkj  e-"  ^  - 1)  (Xr<7'  (3-1) 

t  /  '+/-r 

g-c(r-Tt)^g(c+7®/fc)(/7’+7)  _g(tM-y(B*)/7’^  (3-la) 

it  /  ^+J— 
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B,  when  /  =  -1 


B-1  /  =  -l,  0</<7’ 

Integrate  the  second  term  of  (t)  in  (1-3) : 
MU 

X  bk-\ 

k^-MIl 

M) 

t 

=  JP'Lbk-\  I 

*  ■  -T»y 

t 

=  JP'Lbt.i  f  e<"'>*>*-dX 

*  ’  -r«. 


-00 


=  •«/?  S  *t,-i  e""  -  e<">>tK-J'«»)) 

it  ’  1+7“ 

0<t<Xk 

(3-2) 

it  1+7“ 

0<t<Xk 

(3-2a) 

^■2  /  =  -l,  Tk<t<T 

Mf2  " 

VP  E  6t,-iJ  A(/-x,)e^‘(>-'*)P7<>.-(-r+x*))«/>. 

ife^-M2  4o 

*H0 

«ri_ 

=  i/PE6t.-i  f 

*  -T»y 

=  JP'I.bk-\  e^' 

k  ’  1+7“ 

Xk<t<T 

(3-3) 

=  Jp1L  (1- 

it  1+7“ 

Xk<t<T 

(3-3a) 
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c  wbeiLlzl 


C~L  1^0,  0<t<lk 

PO 

V?  Z  *it.0  J  Kt  -  Pt{^  -  It)  =  0  0<t<Xk  (3-4) 


C-2.  l  =  0,tt<t<T 

Integrate  the  third  term  of  (t)  in  (1-3) : 


oo 

JPUbifi  f 

k  io 

t 

= VF  s  6*0  f  ^y®***^*' 

*  ik 

^■fP'L  6*.o  e^'  ^ 

*  i+y-r 

=  ‘/p'L^ 

*  1+y-e- 


Xk<t<T 

Xk<t<T 


(3-5) 

(3-5a) 


/>.  iS^ 


uUiliU 


UoL-S^ctiti 


Define:  1l{t,.t,)  (SO, 

where 

is  the  unit  step  iunction. 


t\<ti 
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then  K  (tj,  t2)  \.Q,  otherwise 

is  a  unit  pulse  between  tj  and 

hence  for  0<t<T  5^^ (t)  can  be  written  as  : 


^ACI  (^) 

(3-6) 


^  Mil  -1 

■JP  T,  ^  Jt(0, 7)6*,/  e-“  i )  (3-1) 

As^/2  h^L  '  ^+J-r 

JM) 


Mil 


+  Jp  E  Jl(0,x*)6*,-1  (3.2) 

k^n 


hfO 


^  Mil 

+  \jp  Z  »t(T*,7)6*-i  (3-3) 

k^-MTl  ’ 

ktO 


Mil 

*jp  E  n(T*,7)6*.oe-"f^(e<‘^7«..V_e(«-;«>.n.)i  (3.5) 


Mm 

kM 


zl-^ 

(3-1) 

Zl‘-^ 

(3-2) 

zL-i 

(3-3) 

^Xto 

(3-5) 
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IV. 


WQMKASE  ANALYSIS 


A.  Worst-case  for  S(t). 

B.  Worst-case  for  ^  |  S(t)  dt  -  output  of  integrator. 

C  Special  worst-case  analysis. 

D.  Special  worst-case  at  the  output  of  the  integrator. 


A.  WQTstrm^  for  Sit) 

A-1  Worst-case,  ISI : 

From  (2-2),  we  have  : 

=  JP  b.  e-*‘(e‘^  - 1)  E  e‘^ 

1=3-00 


=  JPb.e-^(e‘^- 1)  £  (e-^O* 

h  =  -i 

|e-"|<l 

=  JPb^-“ 

0<t<T 

(4-1) 
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A=2  Worst-case  ACI : 


^k.-i  ^k.0  b,  L  ^ 


r-  ^2  ,  -2 

^cAt)  =JPb  X  { Ji(0,  T^pif  (e(«>«>»)r  _  i )[  J  j 

h^/2 


+Jpb  X  {^c(0,T;t)— e-ic+jm)T)  j 
it  *+i— 

+Jpb  X  {^c(T/fc,  7)-~- 

it 

+  jc(x*,  7)— ^s-  } 

i+/-r 


=  ^fPb  X  «(0,  Th-^sr 

k  *+/-r 

+V?Z^  X  ?c(0,  XitHsr 

*  *-y-r 

X  ^c(Xit,  ly-^  (e/o»t(/-t*)  _  e-c(/-c*)  ^-<ct?to*)r) 

it  »+7-r 

=  X  ^c(0, 7)—^  (g-c(/-^*)  g-(o+/0)t)7’^ 

it  *+y-r 

+/p  6  2  >1(0,  TV-tt  (e*'®^'^*)  -  g-OH?®*)?-) 

t  l+Z-r 

=jpbj:mT)^ 


(4-2) 
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I 

B,  Worst~case  for  ^  J  |5(0P  dt  -  output  ofmte2rator 


Output  of  photodetector  :  ^  \  S(t)\^ 


\S(t)\^  =  \s,(t)  +  SJt)  +  (t)  \ 


=  +  \SJt)\^  +  ftc;  (t)  \'  +  (t)  S]si(t)} 


+ 


2Re{S*B(t)  S^^/t) }  +  2Re  { S]si(t)  S,Jt)  }  (4-3) 


B~1  From  (2-1) 

IS7  (t)  P  =  \SM  =  A  0<t<Z  (4-4) 

T 

/  l5e(/)l  ^dt  =  PreJ  „[1  -  +  2^1  -  e’^0]  (4-5) 


B-2  From  (4-1) 


|5SrO  r  =  Pbf 


0<t<T, 


(4-6) 


(4-7) 
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B-3  From  (4-2) 


MU  M2 


k=-M/2  m=-M2  ( H^) 


M  mM 


=  />l6|2I2 


TH^  (i+>-r)(i-y^) 


0<t<T, 


}  ISS,{0l^<*  =  Pl6P  2  2 

0  km  0+7-r)(l-y-r)  J 


P7lA|2 


=  { 


£  S  f'*- r— .-  1). 


(0 


k  l-K-r)*  ©*  =  ®m 


m 


(4-8) 


(4-9) 


(4-10) 
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2Re  {5«(0 (0>  =  2P*o.o6-(e-" - Q<t<T 


(4-11) 


PTbtifib- 


(1 


(4-12) 


2Re  {5^(0  5^(0}  =  IPb^fiRe  {b  (1  -  e^OS  ^  } 

t  l+y-r 


-  2Pio.o^e  {bl.^ 
k  *-y— 


j  Q<t<T 


n  k  c  y  c  e 


2Re  {S^(0  5^^/)>  =  2Pb-Re  {*  S  ^  q  < ,  <  r 


(4-15) 


j 2Re  {Sr^iO  5;^(0}rfl  =  ^  ^  {*  S  ^  (1  -  (4-16) 
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C  Special  worst~case  analysis 

(dk  =  2nkI/T,  I ;  positive  integer 

T^  =  0 

bo,i  =  b-  €  {0,1} 
bkj  =  bk,o  =  b  E  {0,  eJ^} 

Re{b}E  {0,  cos(t» } 

l*P  e  {0,1  } 

T 

J  I52i(/)I  ^dt  =  PT  bio  [1  -  ^(J  -  +  57(1  -  (4-5) 

0 

Qi  I  l5®(0l  Ut  =  ^  (1  -  (4-17) 

0 


Q3  in  (4-9) 


\\^cAt)\%dt  =  Pnb\^Y.-^ 

n  k  1-K-#)^ 


(4-18) 


C-#  ^ 


/ IRe  {SB(f)S^ (0  (1  - 


0 


(4-12) 


Qj  jzRe  {si(.t)s:ij(t)  }s,dt = (^-ct  _  ^  (4.19) 

n  i  ln~} 


C-6 


2Re  {S^(.t)S:Ut)  )scdt  = 


cT 


k 


(4-20) 
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0 

=  ^  [all  terms  in  section  C]. 

+  [li|2+^i^^^{A}(6.-Ao,o)]Z-;irr}  («l) 

k  cT  ' 
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Now,  we  are  to  encounter  the  most  complicated  situation,  for  clanty,  we 
first  list  the  construction  talbe,  for  which  we  are  to  follow 


S(t) 


2Re{Sl{t)S^a^t)}  2/?e{5a(/)5^(/)}  l5^cy(/)l'  l^/s/Wl' 

I 

1 

C-6 _ ^C-5 _ C-A _ C-3 _ C-1 _ C-1 

I  I  I  III 

2«e{z;‘z';}  2Re{ZlZ^,}  2Re{Zlzh  IZf.!'  IzJI' 


C-3-3 _ ^C-3-2 _ ^C-3-l 

I  I  I 

labteS-Z  2Re{Zl^_,zZ}  \zlj^  \zlj^ 

We  first  calculate  each  tenn  in  Table  (5-2),  then  integrate  it.  For  each  term 
belongs  to  we  also  discuss  two  cases:  (Oii[  =  Q)m  and 
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Since  the  computation  is  very  complicated,  here  we  just  write  down 
the  results  for  each  term,  and  omit  the  intermediate  steps. 


A. 


(4-5) 


B.  From  (2-2) 


-j 


l5/a(0l^=Pe'^'[  S 

/=-Z,o 


0</<7’ 


(5-1) 


^  -1 

J  \Sisi(t)\ ^dt  =  §(1  -  2  (5-2) 

0  i=>-Z,o 

C  Calcualte  -  refer  to  (3-6) 

Q1 


IZJP  =Pe-2-^n:(0, 7)1  ^  -  \)\^  (5-3) 

k  i  *+7— 

^  -2 

lz5'IV/  =  ^(l-  (5-4) 

0  k  I  '  1+7— 
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Q2  '^S  =  min  (t^  T„) 


iz;‘|2  = 


k  m  (l+y— )(l-y— ) 


lgic-ja„)l . 


;fc  m  (‘-‘■y-rKl-y— ) 


j-{C-J<Ok)tgiC-J<Om K-TVtm)  ^-2ct^(cij(iik)(-T-i^k)g(c-J<Om )(-7Vc„ ] 


lzj*lv<=  J  : 


(a)  k^m,  (Ojt^^O) 


ffl 


/T  V  V  L  L* 


[•sk;  -  » 

J\  ^  c  ) 


H — !s;r  - 1) 

H L  (g-(o-y<D*)x,  _  g(c-7®«K-7’+t„) 

_i  -  1)  ^(<>t^WfcK“r4^k)g(<>-yto„)(-r-«„)j 


^(c+yWjtK-T'+Xi)^ 

V  g(c+yo)t)(-7’+tjt) 

']  (5-S) 


(5-6a) 
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(b)  k  =  m,Gik  =  (Om 


%  2  7^  [CT*  +  _  i)e(ct/<..)(-7-«,) 

k  c  )  c 


+_J_-  _  l)g(c->i)(-7+Tfc) 

1-7^ 


_i(g-2«i  _  g(c+>t)(-r-fTt)  g(,c-j<i}k)(-T+Xk) 


+2i?e{-!^(e-<‘^*)^-  ^(c+;tOi)(-r+Tt)^jj 

^■y~ 


(5-6b) 


C-3  3  terms  for  I 

C-3-1  Tg  =  max(Tjt,Tm) 


i'  =  />  2  2  it(x«.  7) 


it  « 


(i+;^Xi-/^) 


Jg(c+/(Dt)t*Q  _  g-(‘^-7®m)7’^j 


(5-7) 


(a)k^m,  G}k^(o 


m 


_PTYY  /--2cx 

cT  ^  ^  0 


it  m 


^  _  ^-2cr\ 

^  (l+7?Xl-7^) 


.)r)j 


(5-8a) 
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ib)k-m,  0)^  =  0)„ 


_  ^  y  IV-il^ 
-  2cr  f 


(1  -  e2(-c7’-fcT,))[j  ^  ^-2cr  _  2e-^^Re{ej'^^'^}] 


C-3-2 


ifc  m  (1+y— )(i-y— ) 

JgjCat-iBm)/  _  g-(c+i<amy^  ^(c+7<B*)tt  _  ^-{c-jiiy^)t  ^{c-jam)Zm 

^-2cf  g(c+;fi»  jt)Xk  g  (c-;<o  „)x«  j 


/lZ?;jVr=JlZ?;.Prf/: 


{a)k^m,  (Hk^dim 
_  fT’  V  ^r*  L  L* 

-  ,J.  Z-  Z-  »*,ofi„,o 


[ 


A:  w 
I 


gjc^kH 

l+J^ 


(gj(0)ik-Om)7’  _  gj(<ak-<am)x*) 
(g-(c+;®m)7’  _  g-(c+;«)w)T,J 


g(o-/»«)tm 


_i(g-2cr  _  g-2CTg  J  ^(c+/<Bk)Tk^(C-7©«)T«  J 

ib)k  =  m,(Ok  =  (iym 
=  f  S  ^  .[cr-  CT*  + 1(1  - 

it  c  ^ 


+2Re{^^ 

i+j-f 


(g  -(c+7®  k)T  _  g  -(c+;a  t>rk)  J  ] 


(5-8b) 


(5-9) 


(5-lOa) 


(5-lOb) 
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C-3-3 


2Re{Zl,.,Zr,, }  =  IPReil.  2 

k  n 


[g-(C+/“)mX  _g-(C+/‘0t)7'^ 

Q(c-j(am)^m  g(C+/©t)Ti^|  _  g-(^+y©/t)7’^JJ 


■2^{Z?,..Zj;„}rf/=  j2ite{Z?;..<o}<^'  ,  for  i;4mand/t  = 


m 


=  S  bk-xb*^,  e^:::±  .(i  e 

km'  '"’°(i+y^KH¥)  ^ 

[ — !^(g-(c+y<o»,)Tg  _g-(c+y<o*)7’j 
-i(g-2ctg  _  g-2cr)  J 


,(c+y©*)T* 


(5-11) 


(5-12) 
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£d 


2Re{Zlzl‘  }=2P/ie{'Li:  2  7r(0. 

k  m  l=-L  (1+7— )(1-;—) 

[g-(<^+;“<iiVg(‘^7“*)(^7’+T*)(g(«^+7W*)7’  _ 

— ^-2«  ^(c->«K-7Vc„)^(c+;(i)i)(/7’+Tt)^^(c+;(Ok)r  _  |jjj 


-2 


g-Rmi-X  1.-10  mtm) 


*  ^  •  IW 

,  'iReiZ^Z'^  }dt  =  /  2iJe{Zjz;‘  }rf/ 


for  k?^m  and  k=m 


it  III  /=— i  V*v  c  A*“J  c  ) 


[ 


1+/ 


e(c~7<amX-7'+t«)jj 


(5-13) 


(5-14) 
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C-5 


-2 


2Re{Zlzr, }  =  2PRe{Y>  S  2  7t(T„,7) 

k  m  I 


+**./*;o 

-*w6«,o 


j2^{zfz?:}rfr=  /2^{zjz?:}* 


tM 


=  Z  Z 


-2 


it  m  l^L  (l+2-rKl-y=?^) 


[  (e“2cxm  _g-2c7’)(l  -Q-ic-jiam)T^{c-j&m'nm 


-  (e’^"  - 


l)^(c+;WiK/7'+TA) 


(5-15) 


(5-16) 
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C-6 

2Re{Zl‘‘Zj*}  The  product  is  only  survived  when  Tm  <t<Tk 


-bk-ibl,o 

+b]t.ib*„Q 

(5-17) 
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\2Jie{Zl'zr,}dt=  f  2Re{ZoZr,)dt 


jm.-i ^g-2cTm  _g-(c-y©m)7’y(c+y(B*X-^+Tt)g(c-yo«)Tg 


I  _g-2CTit)g(0^/0)tX-y+t*)g(g-ytom)TwJ  J 


(5-18) 


a  J  2ife{5i,(0  S'^t))  dt  =  \ 258(0  5’jsf(0  dt 


=  2/’Jo,o  2  ^0^  (g('+‘)^J' - e»^)  J (g-^'  -e  ^dt 


i^Lo 


=  ^  ^o,o(l  -  2  ^0,1 

i=—Lo 


(5-19) 
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£  2«e{5;(0  5^a(0} 


=  2PRe{l£  Jt(0,  T)(e-‘'  ~ 

k  I  1+;— 

+Z  7t(0,  Tt)6o  obt  -1  -  (e-^'  -  e-2<^')e(<^»‘:K-7'+xi)] 

it  ’  ’  1+7-r 

+Z  it(Tt.  7)6o.o6*,-if^(e-"  -  e-'")(I  - 

^  ^+yT 


+2  n(Tk,  T)bofibk,o 

^  •+/  e 


l^^kt  _  g-{c-jis^k)t  _  _  ^-2c/)^(c+;(0i.)ti  j  j  ^5_20) 


I 

J  2/fe{5J(0  5’^a(0}‘* 


—2 

=  ^/je{E  2  ^  1)  g(c+/(0iX^7Vtt) 

it  / 


1+y? 


+2  6o,o6*,-.f^  [(^)  -  {^=^)  - 

k  i+J-T  ;-r  1-y-r 

eWt  ^2_  ^-<ct;a)*)7’y(c+ya)t)t*j 


+2  bdfibk-i^^  [( 

it 


g-i^k^k  ,e-<‘^k>»k-e-<‘^^k)^ 


+2io.o*w^[(^=^)-( 


i+;-r 


•* 

J-r 


l-J^ 


) 


2«-«k-2e-«r_g-2ct*^_2cT 


(5-21) 
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2Re{Sisi(t)  SAa(t)} 

i--Lo 

ReiZ  f.  rt(0,  T)bne-^' 

k  l=-L  ’  1+7-r 

+S  7C(0,  tO  ^-t-1  -  e-2"e('^+>*)(-7’+Tt)) 

it  ’  l+y-# 

+2  3t(T/t,  7)  g'2c/(j  _g-<c+>t)7')g(c+ycoi)t, 

it  ’  1+/— 

+2  Jt(T*,  7)  it,0  ^  («-«->*)'  } 

*  >+J-  (5-22) 


;; 

2/?e  {*S'Jsy(0  5'^c/(0  } 


=  ^  ^  *0^  -  e“^)7?e{zf  _  ])e(«».V7'gC. 

i=-Lo  k  I  2il+j-^) 

+  2  [(^^7^)  -  e<^‘V  J 

k  ^+J—  H-r  ^ 

+  2  — 
k  2(l+y— ) 


+  2  [-TT  [e""‘  -  e->>»^»e-<‘^>>‘)^)  -  i(e'"‘  -  1  } 

it  l+/-r  l-;-r  J  ^ 


(5-23) 


VL  EXACT  ANALYSIS  -  SPECIAL  CASE 

Now,  we  put  two  constraints  into  all  tenns  of  exact  analysis  to  form  a 
special  case  ; 


1.  ©*  = 


Inkl 


I :  positive  integer. 


2.  =  T  for  all  adjacent  channel  k. 


For  the  purpose  of  normalization  for  future  use,  we  set  a  =  ,  or  =  x  = 


so  that  (D/tXyt  =  cojtx  =  ^ar=  27tA:/a 


Still,  we  follow  Table  5-2. 


A. 


B. 


I  |5b(0P<*  =  PTbloW  -  ] 


0 

T 


(6-1) 


-I 


f  l5w(/)l  f 

0  /=-Z,o 


(6-2) 


C.  Special  case  for  ACI 


C-1 


-•7 


k  l=-L 


=  £2^(g2ou:r_g2(a-l)c7-)(gCr_  |)2  | 

^  k  I 


[2 

cr+y2Jt/k/ 


(6-3) 


aT 
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01 

T  X 

JlZJIV/  =  Jlz;Pd/: 

0  0 


(2i)k^m,Gik  ^  COw 

=  (.cT)^PT'Z'Z 


bk.-\  *m.-l 


[ 


k  m  (cT+j2itkl)(cT-j2nntl) 
j2nl(k-m) 

g-c7_g(a-l)cr-y2iwi/a 

cT+j2nml 

g-<T_,g{(».-l)cT-j2i(kIa 

cT-jlitkl 


p2(pL-l)cT_g-2cT  _ 

+- - - -  1 

^  2cT  ■* 


(b)^  =  /il,(Dit  =  ©/f» 


^{cTfPTY. 


I6 


VI I 


(c7)M2wU) 


(p2(»-l)e7'_^-2cT) 

2  2cr 


+2/?e{ 


cT+jlnkl 


}] 


(6-4a) 


(6-4b) 
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C-i  3  terms  for  IZ^P 


C-3-1  j  =  ]  \zl,,Pcll 


(a)  m,(Ok^  (ii 


m 


=  ^cTPT(l  -  S 


m,-! 


k  m 


(cT+j2iUd)(cT-j2nmJ) 


(b)  k  =  m,(Ok  =  (n 


m 


=  \cTPT{\  -  e-^nx  - 


\\Zla\^dt  =  \\Z%\^dt 

0  t 

(a)  k*m,<Sik*fS>m 


=  {cTf  PT'E'Z 


bkfibl,Q 


k  m 


(cT+j2itkI)(cT-j2nml) 


[ 


(.e 


-j2'alcu(ik-m)_ 


1) 


jlid{k-m) 


+ 


^g(a-l)c7'-ty2TO«/a_  j  ^ 

cT-kfliml 

^^{a.-\)cT-j2TtkIa._Y^ 

cT-flnkl 


(l-e2(a-l)c7’) 


2cT 


(6-5a) 


(6-5b) 


(6-6a) 
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(b)  k  =  m,  cO/t  =  co 


m 


=  (c7)2prE 

k 


(c7)2+<27Dfc/)2  2cr 

.-.r^  ,  ^  T 

+2^^  1 


C-3-3 


W/=j2/Je  {Zj;.,Zf„ 


(a)  k^m,(0ic*(0m 


=  KcTfPTRe  {E  S 

jk  m 


{cT+j2vkr)(cT~J2izml)  ' 


^j_^(a:-l)cr-y2iim/a^ 

^  cT+j2nmI 


+ 


g2(ct-l)cT_i 

2cr 


]} 


0>)  ^=/W,  0)t  =  ®m 


=  2{cTfPTRe(Z 

k 


(c7)^+{2wk0^ 


^(l_e(a-l)cr-V2i«Wa)  g2(o-l)cr_j 

L  cT+j2itkI  '*'  2cr 


]} 


(6-6b) 


(6-7a) 


(6-7b) 
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/  2Re  {  Zlr^  )dt  =  ]‘2.Re{  ZlZ]'  )dt 


k^m,  Gik^a^m 


=  2(c7)^TO{SI2^|^|^ 


k  m  I 


_gacT-^miIai  --cr_-(2o-l)c7’  _  ^ 
2CT  ]} 


(6-8a) 


(b)  k  =  m,  Gik-dim 

=  2ic7)^PTRe{rf-^^^^p 

^gacT*fliAla_^  ^  g-cT_g(2a-l)cT  ^  ^ 
I-  cT+j2ickI  57r  J} 


C-5 


i  I 

J  2Re  {  Zlzr  }dt  j  2Re  {  ZlZ^  }dt 


(6-8b) 


=  2{cTf-PTRe  {222 


g{MyT_glcT 

icT+j2nkI)icT~j2nmI)  >-  2cT 


bkjblo  (e2«*-i)-^-l) 

cT+j2nmI  2cT  ^  ^ 


(6-9) 


C-6 


/  T 

J2;je  {Z;zrM/=  j2iJe{Z5Zr}*  =  0 


(6-10) 
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I 

j2Re{SB(t)S;s,(t))dt 


PTbo,o 

cT 


i=-Lo 


(6-11) 


1 

\  2Re{  5J(0  SACiit)  }dt 


2PTRe  2(c7’+/2iti/) 


bQfibk,-icT  (e~^T_^-flidda>^  (j.g-^xcTjjgCo-iyr 


cT+jliUd  >-  fhdd 


cT-flnkl 


6o.o6jt.-i 

2(cT+j2itkI) 

bofibkjocT  ,.(g-y2««^°-i)  (g-c7--y2icfc/a_g-ac7-)  (2-2e(°~^>^^-e~^^-f^<°~^^^) 


cT+flidd  *•  72itjk/ 


cT-flnkl 


1 

F-  j2Re{S'jsAt)SAa(f))dt 


=  2PT[  2  -  e“0  ] 

i=-Lo 
V  *jy 

k  I  2icT+j2nkr) 


yi  6*,_j<r7’  (g(«-i)cr_^-(«+iX7’) 

cT+jlnkl  >■  cT-j2nkI  2cr  -I 


(6-12) 


6jt.-i 

2(cT+j2nkr) 

bkfiCT  (g-«*^7'-e<«-^)^^). 

cT+j2nkJ  ^  cT-j2nkI  2cT  ^ 


(6-13) 


64 


APPENDIX  B 


Programs 


%  exact  analysis  with  optimal  threshold  signal  out  of  the  integrator,  the  signal  X  contains  desired 
signal  and  ACI  &  ISI  and  postdetect  noise  ,The  formula  for  the  BER  is  : 


^^here  p(b)  =  kiif‘i-T)^ 


2%  T 

[J.J  •^un  't-MQ  ■tAf/2) 

MO  MO 


-M2 . . . d^Mti  dx -MU -.diMn 


In  T 


+  /.../  / .../ . . .d^Mii  dx-Mn—diMn] 

MO  MO  i/"o^ 


M  is  the  number  of  adjacent  channels.  Here  we  assume  all  channels  have  the  same  phase  <j>  and 

time  delay  X  i.e,  =  <|),  Tit  =  X  hence  for  our  model  the  power  of  (1/27C)  and  (1/T)  are  fixed 
to  1,  therefore  we  only  have  two  integrals  and  two  arguments. 

M=4;  k=[-M/2:-l  l;M/2]; 
m=k; 

%  produce  the  controlled  matrix  b  to  control  64  different  bit  patterns 
ml=[  zeros(l,32)  ones(l,32)  ]; 
m2=[  zeros(l,16)  ones(l,16)  zeros(l,16)  ones(l,16)]; 
m3=[]; 
for  i=l:4 

m3=[m3  [zeros(l,8)  ones(l,8)]]; 
end 
m4=n; 
for  i*l  :8 

m4=[ni4  [zeros(l,4)  ones(l,4)]]; 
end 

for  i-l:16 

m5=[m5  [zeros(l,2)  ones(l,2)]]; 
end 
m6=[]; 
fori*l;32 

m6=[m6  zeros(l,l)  ones(l,l)]; 
end 
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b=(m  1  ;m2;m3;m4;m5;m6]’; 


%  signal  to  noise  ratio  range  in  dB 
RPTN_DB=[10;0.5;25];  %RPTN_dB 
ppp=10/X0. 1  ♦RPTN_DB); 
lenl=length(ppp); 

%  solalph  function  is  provided  by  Professor  Randy  L.  Borchardt 


n=10; 

[bpx,w6c]=grule(n);  %bpx=bpy  ,w6{=wfy 

%  single  channel 
BER0=0. 5*erfc(ppp/8'^.  5); 

pp»[];  thiesh=[];  %thresh  is  not  normalized  threshold 

for  CT=[5  10  IS  20]  %  cT  is  Fabry-Perot  filter  parameter 
ifCT=5 
qqq=[2:12  17  20]; 
dsdfCT=10 
qq(rt3;12  17  20]; 
dseifCT==15 
qqq=[4:12  17  20]; 
dseifCT**20 
qqq=[5:12  17  20]; 
end 

pe*=0;  threshl*[]; 

for  I=qqq 

BER=[];  thresh2=Q; 
forRPTN=ppp 

tq>=linspa(^0, 1 82, 1 1 );  %  q)pro:dinated  thresholds 

x3=0;  %  firet  few  loops  to  find  out  the  threshold  which  make  the 

%  BER  minimun  dont  care  about  the  scale 

for  i=l  :64  x3=x3+... 

solalph('xx00',0,2*pi,2,bpx,w6c,0, 1 ,2,bpx,  wfic,ap,b,CT,i,I,k,m,RPTN)+. . . 
solalphCxxir,0,2*pi,2,bpx,wft,0,l,2,bpx,wfic,ap.b,CT,i,I,k,m,RPTN); 

Old 

[val,ind]=inin(x3); 

lef=iq}(ind)-16.6; 
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if  lef<0  ,lef=0,  end  %  to  aviod  the  threshold  go  beyond  the  negative  side 


ap=linspace(lef,ap(ind)+ 16.6,1 1); 
%  two  more  time  to  find  alpha 
x3=0; 


fori=l:64  x3=x3+.. 

solalph('xx00',0.2*pi,2,bpx,wfx,0, 1 ,2,bpx,wfic,ap,b,CT,i,I,k,m,RPTN)4-. . . 
solalph('xxl  r,0,2*pi,2,bpx,w6c,0,  l,2,bpx,w6c,ap,b,CT,i,I,k,m,RPTN); 
end 

[val,ind]=min(x3); 

lefi=ap(ind)-3.1; 

if  lef<0  ,ld==0;  end 

ap»iinspace(lef,ap(ind)+3 .1,11); 

%  three  more  time  to  find  alpha 
x3=0; 


for  i=l  ;64 
x3“x3+... 

solalph('xx00',0,2*pi,2,bpx,w6c,0, 1 ,2,bpx,  w6c,ap,b,CT,i,I,k,in,RPTN)+. . . 
solalphCxxl  r,0,2*pi,2,bpx,w6c,0,  l,2,bpx,wfit,i^),b,CT,i,I,k,m,RPTN); 
end 

[val,ind]=mm(x3); 

lef=ap(ind)-0.57; 

if  lef^O  ,lef=0;  end 

ap=linspace(letap(ind)-K).  57,8), 

%  four  more  time  to  find  alpha 
x3*0; 


for  i=  1:64 
x3=x3+... 

solalph('xx00',0,2*pi,2,bpx,Av6c,0, 1 ,2,bpx,w&,ap,b,CT,i,I,k,m,RPTN)+. . . 
solalphCxxl  l',0,2*pi,2,bpx,wfic,0,l,2,bpx,wfic,^,b,CT,i,I,k,m,RPTN); 

end 

[val,ind>’miin(x3); 

ap=‘ap(ind); 
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%  after  find  optimal  alpha  use  double  integration 

^3=0; 

for  i=l;64 

qq3=qq3+(  dbgquadm('xx00',0,2  •pi,2,bpx,  w6t,0, 1,2,... 
bpx,  w6c,ap,b,CT,i,I,k,m,RPTN)+. 
dbgquadm('xx  1 1  ’,0,2*pi,2.bpx,  wfic,0, 1,2,... 
bpx,wfic,ap,b,CT,i,I,k,m,RPTN)  )/(256*pi); 
end 

BER*[BER  qq3];  thresh2=[thresh2  ap]; 

if  qq3<l  0^- 1 5)  %for  save  time  only  interesting  in  1 CK^- 1 5) 

ii=find(ppp=RPTN); 

BER(ii+l:lenl)=5*l(KX-116)*ones(  l,length(ii+l:lenl) ); 
thresh2(ii+l  :lenl)=(ap+2)*ones(l,length(ii+l  ;lenl)  ); 
break 
end 


pe=[pe;BER];  threshl=[threshl;thresh2); 
end 

pp=[pp;pe];  thresh=[thresh;threshl]; 
end 


time2=toc; 


%  worse  case  form  appendix  setting  optimal  threshold  equal  (x0+xl)/2 

I=linspace(0,6, 101);  %  set  I  value  in  linspceinteger  to  find  minimum  1 

RPTN_DB=®0:0.2;20;  %  x-axis  dB  range 

RPTN=10/t0. 1  •RPTN_DB);  %  change  to  ratio 
%  single  channel 
BER0=0.5*erfc(RPTN/8^.  5); 

%find  optimal  alpha 
pp=[]; 

for  CT=[5  10  15  20] 
x3=0; 

fork=l:0.5*M 

x3=x3+2./(l+(2*pi*k*I/CT)/^2); 

end 

xO=(l-exp(-2*CT))/(2*CT)+x3*(l+2*(l-exp{-CT))/CT); 

xl=l-2*(l-^(-CT))/CT+(l-exp(-2*CT))/(2*CT); 

[val,ind]=min(abs(xO-xl)); 

I=ceU(I(ind)); 

I=(I(ind));  %  our  minimum  I  value 

%  set  different  I  and  to  find  the  BER 
ss=I+l:I+10; 

BER=[]; 
for  I=ss 
x3=0; 

for  k=l;0.5*M 

x3=x3+2/(l+(2*pi*k*I/CT)^2); 

end 

xO=(l-exp(-2*CT))/(2*CT)+x3*(l+2*(l-exp(-CT))/CT); 
xl=l-2*(l-exp(-CT))/CT-Kl-exp(-2*CT))/(2*CT); 
BER=[BER;0.5*erfc(  RPTN*(xl-x0)/8'X).5  )  ]; 
end 

pp=[pp3ER]; 


semilogy(RPTN_DB,BERO.'--',RPTN_DB,pp(;,:)) 
axis([10  19  l(yx-15)  1]) 
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